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Abstract 

We study how the curvature of spacetime, in conjunction with solar radiation pressure (SRP), 
affects the bound orbital motion of solar sails. While neither the curvature of spacetime nor 
the SRP alter the form of Kepler's third law by themselves, their simultaneous effects lead 
to deviations from this law. We also study deviations from Keplerian motion due to frame 
dragging, the gravitational multipole moments of the sun, a possible net electric charge on 
the sun, and a positive cosmo logical constant. The presence of the SRP tends to increase 
these deviations by several orders of magnitude, possibly rendering some of them detectable. 
As for non-circular bound orbits, the SRP dampens the rate at which the perihelion is 
shifted due to curved spacetime, while the perihelion shift due to the oblateness of the sun 
is increased. With regards to the Lense-Thirring effect, the SRP increases the angle of 
precession of polar orbits during one orbital period, although the precession frequency is not 
actually altered. We also consider non-Keplerian orbits, which lie outside of the plane of 
the sun. In particular, we investigate how the pitch angle of the solar sail is affected by the 
partial absorption of light by the sail, general relativistic effects, and the oblateness of the 
sun. Non-Keplerian orbits exhibit an analog of the Lense-Thirring effect, in that the orbital 
plane precesses around the sun. A near-solar mission for observations of these effects could 
provide an interesting confirmation of these phenomena. 
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1 Introduction 



In the last decade, the observation and analysis of satellite motion has provided an abun- 
dance of data with which to test basic physical principles. Examples include the Pioneer 
anomaly, which is an unexplained acceleration of the Pioneer 10 and 11 spacecraft on es- 
cape trajectories from the outer solar system [1,2], and the flyby anomaly, for which the 
velocities of the Galileo, NEAR and Cassini spacecraft are different from what is expected 
after Earth flybys [3,4]. In fact, the difficulties of explaining these anomalies within the 
framework of standard physics became a motivation to speculate on the unlikely possibility 
that they originate from new physics. Missions have been proposed [4-6] to further explore 
these anomalies, in order to better understand the laws of fundamental physics as they affect 
dynamics within the solar system. 

One of the most basic laws that describes motion in the solar system is Kepler's third 
law, which can be derived from Newton's law of gravitation and provides a relationship 
between the period T and the orbital radius r of an object orbiting the sun. Namely, T 2 
is proportional to r 3 with the proportionality constant given by 4ir 2 /GM, where G is the 
gravitational constant and M is the mass of the sun. Below, we discuss deviations from 
Keplerian orbits due to phenomena within conventional physics, which can be observed from 
the motion of solar sail propelled (SSP) satellites [7,8]. 

In the Newtonian approximation, the sun is the source of a gravitational force on other 
masses. In the general relativistic framework, in the absence of non-gravitational forces, 
objects follow geodesies on the curved spacetime in the vicinity of the sun. At the same 
time, the sun is also a source of solar electromagnetic radiation, which produces an external 
force on objects via the solar radiation pressure (SRP). It is of particular interest to analyze 
how the trajectory of an object deviates from a geodesic under the action of the force due 
to the solar radiation pressure. We will assume that the backreaction of the radiation on 
spacetime is negligible. Therefore, we can say that objects move in the photo- gravitational 
field of the sun. 

When an SSP satellite undergoes orbital motion within the photo-gravitational field of 
the sun, the orbital period is altered by the presence of the solar sail, since the force from the 
solar radiation pressure affects the dynamics of the orbit. Moreover, due to the continuously- 
acting solar radiation pressure, an SSP satellite is capable of exotic non-Keplerian orbits 
which are impossible for spacecraft which are not equipped with solar sails. The solar sail 
enables one to design trajectories for which the SSP satellite is effectively levitated above 
the sun- namely, the sun lies outside of the orbital plane. For example, non-Keplerian orbits 
can be above the ecliptic plane and parallel to it. 

Satellite orbits been extensively studied within the Newtonian approximation for gravity 
(see [9], for example), and many papers have considered general relativistic effects on orbits. 
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However, general relativistic effects on solar sails, which are propelled by electromagnetic 
radiation pressure, has largely remained an unexplored subject. At first glance, such an 
intersection of topics might not appear necessary. However, we currently have the technology 
to put solar sails in orbit around the sun with heliocentric distances approaching 0.05 AU. 
At such a proximity to the sun, which is considerably closer than the orbit of Mercury, the 
continuous effects of curved spacetime on bound orbits should be considered. 

We recently considered various general relativistic effects on solar sails in circular orbits 
within the plane of the sun [10]. The curvature of spacetime can be separated into a contri- 
bution associated with a static central mass and a portion corresponding to the rotation of 
the sun. For the purpose of comparison, we also considered a number of other effects, such 
as the gravitational multipole moments of the sun, a possible small net electric charge on 
the sun, and a small positive cosmological constant. 

In this paper, we generalize this study to include slightly non-circular orbits, orbits at 
arbitrary polar angle (since spherical symmetry is broken by frame dragging effects), and 
non-Keplerian orbits. We find that the perihelion shift of non-circular orbits which is due 
to curved spacetime in general relativity actually occurs at a slower rate due to the SRP. 
On the other hand, there is also a perihelion shift that occurs in the reverse direction due 
to the oblateness of the sun. This second perihelion shift is also present within Newtonian 
gravity, and we find this effect to be augmented by the SRP. We also consider the Lense- 
Thirring effect, which is the precession of polar orbits due to frame dragging, and which has 
no counterpart in Newtonian gravity. Perhaps surprisingly, the precession frequency is not 
altered by the SRP. We also consider various effects on solar sails undergoing non-Keplerian 
orbits, such as the partial absorption of sunlight by the surface of the sail, spacetime curvature 
and frame dragging, and the oblateness of the sun. Interestingly enough, we find an analog 
of the Lense-Thirring effect for non-Keplerian orbits, for which the orbital plane (which lies 
outside of the plane of the sun) precesses around the sun. 

This paper is organized as follows. In section 2, we consider various effects on bound 
orbits in the plane of the sun, including solar radiation pressure, static curved spacetime, 
frame dragging, the gravitational multipole moments of the sun, a possible small net electric 
charge on the sun, and a small positive cosmological constant. In section 3, we consider 
circular orbits out of the plane of the sun. After a brief review of non-Keplerian orbits for a 
perfectly reflecting solar sail, we consider the effects of light absorption, spacetime curvature 
and frame dragging, and the oblateness of the sun. Conclusions follow in section 4, where we 
summarize our results, and consider a solar-sail propelled satellite with specified structural 
and orbital parameters in order to do a comparative analysis on the how the SRP affects 
various phenomena. 
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2 Bound orbits in the plane of the sun 



The purpose of this section is to point out various sources of deviations from Keplerian 
orbits and discuss how the resulting change in period is enhanced by the SRP to the degree 
in which it may be observed for some cases. The phenomena discussed include the curvature 
of spacetime in the vicinity of the sun, described by the Schwarzschild metric, frame dragging 
due to the rotation of the sun, for which the curved spacetime is described approximately 
by the large- distance limit of the Kerr metric [11], the oblateness of the sun, the effect of a 
possible small net electric charge on the sun, and a small positive cosmo logical constant. This 
list is by no means exhaustive. Some effects that we do not consider are the perturbations 
in motion due to the gravitational fields of the planets, the magnetic field of the sun and the 
solar wind. While the reality is that all effects occur simultaneously, particular effects may 
be isolated by considering a variety of orbits. However, in order to get a fairly accurate idea 
of the relative importance of some effects in the presence of the SRP, we will consider them 
individually for the simple scenario of circular orbits. 



2.1 Solar radiation pressure 

We will be mainly considering an SSP satellite orbiting the sun. According to Maxwell's 
electromagnetic theory, electromagnetic waves carry the energy and linear momentum and 
the radiation pressure P exerted on the surface of satellite or solar sail due to momentum 
transport by photons is given by 

where S is the magnitude of the Poynting vector and the solar luminosity is L s = 3.842 x 10 26 
W. Also, 0.5 < rj < 1, where rj = 0.5 corresponds to the total absorption of photons by the 
surface of the satellite and rj — 1 corresponds to the total reflection of the solar radiation. 
The resulting force on the SSP satellite is F = PA, where A is the area of solar sail directly 
facing the sun. Thus, the acceleration due to this force can be expressed as 

K= w a= A' (2) 

where m is the mass of the SSP satellite. The mass per area a is a key design parameter for 
solar sails [7, 12, 13]. 

We will first consider the effect of the SRP on Keplerian orbits in the Newtonian ap- 
proximation for gravity. The SRP force is repulsive and the gravitational force is attractive. 
For simplicity, we will start by restricting ourselves to the case in which the surface of the 
solar sail is directly facing the sun, so that vectors normal to the sail are directed along the 
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Figure 1 : The gravitational force and solar radiation pressure force acting on a solar sail which is directly 
facing the sun. 

sun-satellite line. Then both forces act along the same line, as shown in Figure HJ Also, 
both forces fall off as 1/r 2 , since r is the heliocentric distance. Therefore, the consideration 
of both forces leads to a modification of the effective mass of the sun in Kepler's third law. 
Namely, the mass of the sun, which is M = 1.99 X 10 30 kg, can be effectively renormalized 
as M = M-k/G. 

The modified Kepler's third law can be expressed as 

T 2 = i!tr 3 , (3) 
GM 

where r is the radius in the case of circular orbits, while for elliptical orbits we replace r by 
the length a of the semi-major axis of the ellipse. Since M < M, the orbital period of an 
SSP satellite will always be longer than that of a conventional satellite for a given orbital 
radius. In the Newtonian approximation, Kepler's third law retains its form in the presence 
of the SRP. Eq. (j3J) can easily be obtained from the usual expression for Keplerian orbits by 
simply replacing the solar mass M with the effective solar mass M. This is just a reflection 
of the fact that both the Newtonian force of gravity and the solar radiation pressure force 
go as the inverse square of the heliocentric distance. One can show that the orbital radius r 
(length of the semi-major axis a for elliptical orbits) can be expressed in terms of the energy 
of the SSP satellite. Thus, the period can be expressed in terms of the energy as well as the 
solar sail parameter k. We can also use to express k in terms of the period and radius as 

4^2 3 

k = GM , (4) 

which is useful in the design of SSP satellites. 

As a first example, we will consider Mercury, whose average orbital radius is r = 5.79 x 
10 10 m, which corresponds to a period of about 87.9686 days. Mercury has a mass of 
3.30 x 10 23 kg and a radius r M = 2.44 x 10 6 m. This yields o = 1.76 x 10 10 kg/m 2 , where 
the effective area is itr 2 M . As a first approximation, we assume that no sunlight is reflected 
by Mercury, so that rj = 0.5. This yields an increase in period on the order of 10~ 7 s which, 
as to be expected for any planet, is negligible. 

We will now consider a conventional satellite orbiting the sun at a distance of r = 1 
AU~ 1.50 x 10 11 m, which corresponds to a period of one year. If the mass of the satellite 
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is 1000 kg and its area is 2 m 2 , then a = 500 kg/m 2 . If we suppose that rj = 0.75, then the 
increase in period due to SRP is about 36 s, which could be observed. 

In the remainder of this section, we will consider an SSP satellite with the following 
specifications: 

r = 0.05 AU « 7.48 x 10 9 m, 

r, = 0.85, a = 0.00131 kg/m 2 . (5) 

If the acceleration due to the SRP is ignored, then the corresponding orbital period would 
be about 4 days. When the SRP is taken into account, the period is about 70 days. 

Since we are considering motion in a plane, we can use polar coordinates (r, 0). We can 
derive the orbital equation from the following expressions for conserved quantities 

1 .o 1 2 ; 2 GM 

L N = r 2 <P, (6) 

where dots denote time derivatives, and En and Ln are the energy and angular momentum 
per unit mass of solar sail, respectively. The subscripts N are used to distinguish these 
conserved quantities from their general relativistic counterparts. Then the orbital equation 
is given by 



\d(j) J \ r r 2 J L\ 



n 



and a solution is 



GM ^2E N + G 2 M 2 /L 



r l = -^ + ^ cos0. (8) 

L N Ln 

2.2 Static curved spacetime 

We will now consider the simultaneous effects of the SRP and curved spacetime in the 
vicinity of the sun. We assume that the backreaction of the electromagnetic radiation on 
the background geometry is negligible so that it acts on the SSP satellite only via the SRP. 
The exterior spacetime of the sun is approximately described by the Schwarzschild metric, 
which is given by 

ds 2 = -fc 2 dt 2 + f~ l dr 2 + r 2 dQ 2 , (9) 
where the metric function / is given by 

2GM , x 

/ = l--5-. 10 

c 2 r 
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Also, dVt 2 refers to the metric of a unit two-sphere S 2 , which can be written in terms of the 
polar angle 9 and azimuthal angle (f) as 



dtt 2 = d6 2 + sin 2 9 d<p 2 . (11) 

Spherical symmetry allows us to orient the coordinate system so that the orbit is confined 
to the equatorial plane at 9 = 7r/2, and thus pe = 0. Since the metric is independent of time 
and the azimuthal direction 0, the corresponding components p t and p^ of the 4-momentum 
are conserved. We define the constants of motion E = —p t /m and L = p^/m, where m is 
the rest mass of the SSP satellite. Thus, 

, mE dr a , m , . 

p' = ^-, p r =m— , p e = 0, p*=-L, (12) 
c z j dr r A 

where r is the proper time. In the absence of the SRP, p 2 = —m 2 c 2 yields 

dr\ 2 E 2 /, L 2 



- 



o 2 + -T /• (13) 



dr J c 2 \ r 2 
Differentiation of ( fl3l) with respect to r gives the radial component of the 4-acceleration 

r d 2 r GM L 2 3GML 2 

« P = *3 + -pr-pr + -s3-- ( 14 ) 



We will now turn on the SRP, so that 



where k is given in (j2J) • Note that even though the coordinate r does not measure the proper 
distance, the surface area of a sphere is still given by 47rr 2 , which means that the magnitude 
of the Poynting vector as well as the acceleration are given by the same expressions as in 
the Newtonian approximation. Equating the expressions for a r given in (Tl4|) and (Tl5|) and 
taking the first integral gives 

dr\ 2 E 2 „ 2GM L 2 



-r = — -c 2 + - ¥ f. (i6) 

dr J c z r r z 



From (fl6l) and the equation in (1121 . we find the orbital equation 

'//V ( E 2 2 2GM L 2 \ r 4 

Note that the SRP reduces the effective mass only in the term which is present for Newtonian 
gravity. In the limit of Newtonian gravity and non-relativistic speeds, E 2 /c 2 — c 2 — > 2En 
and the orbital equation ( ITTj) reduces to (J7J). 
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2.2.1 Circular orbits 



We will first consider orbits which are completely circular, for which 



This yields 



dr d 2 r , . 

Tr = ^ ^ = (18) 



2 _ 4 ( AGM - c 2 r \ c 2 GM 2 c 2 GMr 2 

^\c 2 r-3GMj r ' c 2 r-3GM' { ' 

2 



The above equation for L can be inverted to get 



L 2 + x L 4 - YlG 2 L 2 MMlc 2 

r = ¥ = . (20) 

2GM 



We can write ffTBT) in terms of an effective potential as 
where 



c 2 £ =^-^//(r), (21) 



^/ = ^- — {^Mm)- (22) 

The minimum of V e ff is at the radius given by (|20|) . which therefore corresponds to a stable 
circular orbit. 

Using dt/dr = p t jm and d(p/dr = p^/m, we find for the orbital period T: 



rp2 



2^,3 



4ir 2 r 
GM 



c 2 r - AGM 

l + K 



(c 2 r - 2GMf 



(23) 



If we consider the case of k = in (!23l) . which means that there is no effect from SRP, then we 
obtain the classical Keplerian expression for the orbital period. Therefore, the description 
of the motion of an object in the static exterior spacetime of the sun described by the 
Schwarzschild metric does not alter Keplerian circular orbits. While the solar radiation 
pressure diminishes the effective solar mass, the form of Kepler's third law remains intact. It 
is only through the simultaneous effects of general relativity, via the Schwarzschild metric, 
and solar radiation pressure on the orbital motion of the SSP satellite gives the deviation 
from Kepler's third law by the radial dependent (and thus energy dependent) factor in the 
squared brackets. 

Keeping only the leading correction due to the curvature of spacetime, (f2"3"j) reduces to 



T 2 ~ 



47r 2 r 3 r K 
1 + 



GM 



c 2 r 



(24) 
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For the specification given in (jSJ), we find that this yields an increase in the period of about 
0.6 s. 

Inverting fl24"j) to express the k parameter in terms of the period and radius gives 

GMT 2 - 4vr 2 r 3 
K ~ T 2 + 4vr 2 r 2 /c 2 ' (25) 

2.2.2 Non-circular orbits 

We will now calculate the precession of orbits which are not completely circular. We define 
the coordinate 

1 GM , . 

» = ;--tr- < 26 > 

For a nearly circular orbit, y is small. Thus, neglecting terms in y 3 , the orbital equation (fT 
can be expressed as 

dy\ 2 E 2 /c 2 + G 2 M 2 /L 2 - c 2 2G A MM 3 6G 3 MM 2 ( 6G 2 MM \ 2 



L 2 c 2 L 6 c 2 L 4 tf \ c 2 L 2 

The solution is 

y = yo + ^cos[A;(0-0o)], (28) 
where 0o is an arbitrary constant and 



6G 2 MM 

k = \ 1 — 



yo 



c 2 L 2 
3G' 6 MM 2 
c 2 L 4 k 2 ' 



1 E 2 /c 2 + G 2 M 2 /L 2 - c 2 2GMM 3 6G 3 MM 2 

= jfeV— — z 2 "^ ^ 

Note that the orbit oscillates about y = yo which corresponds to the radius for a circular 
orbit in general relativity, whereas y = corresponds to the radius for a circular orbit in 
Newtonian gravity. The perihelion shift during one complete orbit is 

A(f>= — -2?r. (30) 
k 

For nearly Newtonian orbits, 

. , 6nGM . , 

w , 31 

c z r 

where we have taken L 2 m GMr. Note that the perihelion shift per orbit is not affected 
by the SRP at this approximation. In fact, this is exactly the formula which describes the 
perihelion shift for the orbit of Mercury, which is about 0.43 arcseconds per year. The 
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perihelion shift of a conventional satellite at r = 0.05 AU is substantially larger- about 70 
arcseconds per year. Interestingly enough, this is one scenario for which the SRP would 
dampen the effect, since it increases the period of the orbit, and therefore the time required 
for the perihelion shift A0 to occur. Using the approximate period given by ([3]), we find 
that the rate at which the perihelion is shifted is 



ZGMVGM 

T ~ c 2 r 5/2 • y. 6Z ) 

As an example, for our SSP satellite the perihelion shift is reduced to about 4 arcseconds 
per year. 



2.3 Frame dragging 

The rotation of the sun causes frame dragging, which affects the trajectories of orbiting ob- 
jects. The external spacetime of a slowly rotating body with mass M and angular momentum 
J is described approximately by the large- distance limit of the Kerr metric [11] 

AG T dr 2 
ds i = - f( ? d t> - ^± sin 2 e dtd(f) + zl_ + r 2 d Q 2 , (33) 

err / 

where / is given by ffTUl) . We do not use the full Kerr metric since it does not seem to 
describe the external spacetime of a rotating material body, because it does not smoothly 
fit onto metrics which describe the interior region occupied by physical matter. Since there 
are corrections to this metric in higher-order J we will work up to only linear order in J, 
which suffices for the slowly rotating sun. Note that J > for a prograde orbit with respect 
to the sun, while J < for a retrograde orbit. 

As previously, we define the constant of motion E = —p t /m. The spherical symmetry 
of the Schwarzschild metric enabled us to orient orbits to lie within the equatorial plane. 
However, we no longer have this luxury for the metric (1331) . Thus, there is now a constant of 
motion associated with the component of angular momentum per unit mass that is normal to 
the equatorial plane, and an additional constant of motion associated with the total angular 
momentum per unit mass. Namely, L z sm9 = p^/m, where the z subscript denotes the 
direction that is normal to the equatorial plane, and 

a Tn L 2 . , . 

F r 2 V tan 2 6 

where the Carter integral P is related to the total angular momentum [14,15]. Note that 
P = corresponds to motion in the equatorial plane, for which 6 = n/2 and p e = 0. No 
other orbits can lie within a fixed plane. 
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It will be useful to have the contravariant versions of p t and p^. Inverting the (t, 4>) 
portion of the metric ( 1331) gives 

o« - L a t4> - 2GJ H - 1 (oc\ 

y — Of' {J — An 3 J 9 _ 9-2/1' 

cr/ <r/r d H sin 



which leads to 



. H u mi? 2GmJL z sin # 

p = g tt Pt + 9 t *P4> = -2T m — ' 

^ ^ ^ mL z 2GmJE , . 

p* = ^*p* + 9 t *Pt = - r ^ L o + -zrn-- 36 

r 2 sin V crr^j 
In the absence of the solar radiation pressure, p 2 = —c 2 m 2 yields 



dr J c 2 \ r 2 J c 4 r 

Differentiating this with respect to r gives 



d 2 r GM L 2 + P 3Gc 2 M(L 2 + P) - 6GJEL Z 

+ + -• ( 38 ) 



dr 2 r 2 r 3 c*r 



Turning on the solar radiation pressure, the equation of motion is given by ffTol) . which is 
not altered by linear terms in J. Taking the first integral of a r yields 



/ dr_ \ 2 = &_ c 2 + ^GM_ _ f(L 2 z + P) _ AG J EL z 
\dr J c 2 r r 2 c^r 3 



For a circular orbit, dr/dr = and d 2 r/dr 2 = 0. Thus, from (l37j) and (1381) 



E 



X cJ G 2 Y 



r(c 2 r - 3GM) r 2 V (c 2 r - 3GM) 



Y 3GJ X 



c 2 r - 3GM c I/ r(c 2 r - 3GM) 



3 ' 



(40) 



where 



X = c 4 r 2 -c 2 G{3M + M)r + 4G 2 MM, 

3G M P 

Y = r(GMr - P) + — . (41) 

cr 

We will now consider the special cases of orbits in the equatorial plane of the sun and polar 
orbits separately. 
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2.3.1 Orbits in the equatorial plane of the sun 



For rotation in the equatorial plane (8 = 7r/2 and P 
p^/m, we find the orbital period to be given by 



0), from dt/dr = p l /m and d^jdr 



2^,3 



4ir z r 
GM 



1 + K- 



(c 2 r - AGM) 
c 2 r - 2GM) 2 



1 + 



2v^(i + ggyp) 



2 V Mr 3 / 2 



1 + (c 2 r-2GA/) 2 



(42) 



We have presented T 2 as a product of three factors. The first factor is the same as in the 
Newtonian approximation for gravity in the presence of the SRP, and embodies the fact 
that the SRP effectively renormalizes the solar mass. The second factor shows the deviation 
from Keplerian orbits due to the simultaneous effects of SRP and the curvature of spacetime 
outside of a static central body. The third factor gives the deviation from Keplerian orbits 
from the combined effects of the SRP and frame dragging due to the rotation of the sun. Note 
that the second and third factors in (j4*2i) involve both the solar mass M and the parameter 
k separately, rather than simply in the combination of the renormalized solar mass M. For 
the case of k — 0, which means that there are no SRP effects, (|4*2]) reduces to 



2^,3 



T 



4iT 2 r 
GM 



1 + 



2VGJ 



(43) 



c^VMr 3 / 2 

Thus, in the absence of the SRP there is still a deviation from Keplerian orbits due to frame 
dragging. 

Keeping only the leading contributions due to spacetime curvature in (j4*2]l gives 



2^,3 



T 



Att r 
GM 



1 + 



K 



c 2 r 



1 



2VGJ 



ix/MrV 2 



(44) 



Note that the SRP increases the effect that frame dragging has on the orbital period, since 
it is the renormalized mass M that appears in the last factor. The speed of the outer layer 
of the sun at its equator is v ~ 2000 m/s at the equatorial radius R ~ 7 x 10 8 m, and 
we will assume that the core of the sun rotates with the same angular speed. Therefore, 
J = ^MvR ~ 10 42 kg m 2 /s. Without the effects of the SRP, frame dragging leads to an 
increase (decrease) in the period of only 4 x 10~ 5 s for a prograde (retrograde) orbit. With 
the SRP and the specifications in ([5]), frame dragging leads to a change in the period of 
about 0.01 s. 



2.3.2 Polar orbits 

We will now consider orbital planes which pass through the poles of the sun. For this case, 
L z = and the angle along the orbit is 9. Setting L z = in (HOI) gives 

c 2 GMr 2 

(45) 



c 2 r - 3GM 
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Note that the rotation of the sun does not effect this expression up to linear order in J. 
From the expressions for p t and p 6 in (1361) and ( 134"1) . we find the period of a polar orbit to 
be given by (123]) . That is, up to linear order in J, the period of polar orbits is not affected 
by the effects of frame dragging. 

However, the orbital plane of a non-equatorial orbit will precess in the <fi direction, which 
is a well-known effect of frame dragging called the Lense-Thirring effect. From (1361) . we find 
that the precession frequency for a polar orbit is 

^ = rf0 = 2GJ 
p dt c 2 r 3 

Up to linear order in J, the precession frequency is not altered by the SRP. In addition, 
the angular speed given by (|46|) coincides with that of a particle falling from infinity with 
zero angular momentum. Thus, the angular speed in the direction is not altered by the 
polar orbital motion. Also, the precession frequency does not depend on the location of the 
satellite in the polar direction, although this would not be the case if quadratic terms in J 
were considered. We will use these facts to find the approximate precession frequency for 
non-Keplerian orbits in section 3. 

The angle of precession during one orbital period is given by 

^GJ 

A<p w == , (47) 

c 2 VGMr 3 

which is increased by the SRP due to the corresponding increase in period. For both a 
standard and SRP satellite in a polar orbit at r = 0.05 AU, we find the rate of precession to 
be about 0.03 arcseconds per year. 

2.4 Gravitational multipole moments of the sun 

We will now consider the effect of the mass multipoles of the sun on bound orbits. In 
particular, the dominant higher moment is the quadrupole, which is associated with the 
oblateness of the sun. Note that we are still taking the sun to be a point-like light source. 
Our assumption is that variations in the SRP associated with the shape of the sun are 
negligible compared to the gravitational effects of oblateness. 

2.4.1 Circular orbits 

Working in Newtonian gravity, the external gravitational potential of an oblate spheroid is 
given by [16] 

GM GMJ^ fR\ n , „. 
V = - P„(cos#), 48 

n=2 V 7 



13 



where J n are the multipole mass moments, R is the equatorial radius of the sun and P n are 
the Legendre polynomials. Note that the effective mass in the first term is renormalized by 
the solar radiation pressure, whereas the multipole mass moments are not affected. We will 
consider only the case for which the orbit of the SSP satellite is confined to the equatorial 
plane. Then the acceleration is purely in the radial direction, given by 

GM ( 3MJ 2 R 2 15MJ 4 R 4 \ . . 



r 2 V 2Mr 2 8Mr 4 
where we consider the quadrupole and octupole terms. Upon taking a power series expansion 
in the mass multipole moments, this leads to the relation 



2-3 



T 2 



4ir 2 r 
GM 



3MJ 2 R 2 3MR A ( 6MJ. 
1 + J- ^— 5 J 4 + 



(50) 

2Mr 2 8Mr 4 V M 1 

From a model of the interior structure and of the solar rotation obtained from helioseismic 
inversions, the first few gravitational multipole moments have been theoretically predicted 
to be 

J 2 « 2.2 x HT 7 , J 4 « -4.5 x 10^ 9 , J 6 « -2.8 x 10~ 10 , (51) 

which are in reasonable agreement with previous models [17]. However, the quadrupole mass 
moment of the sun was recently been measured with unprecedented precision to be as much 
as J 2 ~ 9 x 10~ 6 during active phases of the solar cycle [18]. 

First we consider the effects of the sub- leading term in (150]) . which is linear in J 2 . Without 
the SRP, this latest value of J 2 increases the period by about 0.02 s. With the SRP and the 
specifications given in (jSJ) for the SSP satellite, J 2 increases the period by about 105 s. 

We now turn to the next highest terms in the power series expansion, which includes a 
linear term in the octupole mass moment J4 and a quadratic term in the quadrupole mass 
moment J 2 , as can be seen from (1501) . Without the SRP, the linear J 4 term dominates over 
the quadratic J 2 term and increases the period by 10~ 7 s. With the SRP, the linear J4 term 
in (l5Dj) increases the period by 6 x 10~ 7 s. However, due to the SRP, the quadratic J 2 term 
now dominates over the linear J4 term and leads to a decrease in the period by 0.003 s. 
Whatever the case, except for the linear J 2 term, the higher mass multipole terms have less 
of an effect on the period than even frame dragging. 

2.4.2 Non-circular orbits 

From the Newtonian equations for energy and angular momentum, and taking into account 
the mass quadrupole term in the gravitational potential given by (j4"8"j) . we find the orbital 
equation to be 



dr\ 2 E 2 „ 2GM L 2 GMJ 2 R 2 \ r 4 . . 

c2 + — Z5T— 72- ( 52 ) 



- 



14 



Note that this equation has the same form as the general relativistic orbital equation given 
by (117j) . Therefore, orbits which are not completely circular will undergo precession due to 
the oblateness of the sun. Using the coordinate y given by fl26|) . we will consider a nearly 
circular orbit for which y is small. Neglecting the y 3 terms, the orbital equation can be 
written as 

dy\ 2 E 2 /c 2 + G 2 M 2 /L 2 - c 2 G A M 2 M 3 3G 3 M 2 M 2 ( 3G 2 M 2 M , ., 

1 " " y- — TT— + 1 V > ( 53 ) 



d<t>) L 2 L 8 L 6 \ L 4 

where M 2 = MJ 2 R 2 . The solution has the same form as (!28l) . where now the constants are 
given by 



3G 2 M 2 M 

1 + 



2/o 



3G 3 M 2 M 2 



A = ^E^Gp^_G^_3G^M lyo _ k2yi (54) 

The orbit oscillates about y = yo and the perihelion shift during one complete orbit is given 
by ( J301) . For small oblateness, 

Mr 2 K J 

where we have taken L 2 « GMr. The characteristics of the precession due to the oblateness 
of the sun are rather different than the general relativistic precession. For instance, this 
precession occurs in the reverse direction relative to the orbit. For Mercury, we find the 
rate of precession to be about —0.01 arcseconds per year. Note that this is substantially 
less than the precession of Mercury's orbit due to the curvature of spacetime, which is 
about +0.43 arcseconds per year. For a conventional satellite in orbit at r = 0.05 AU, we 
find the precession rate to be about —14 arcseconds per year. For our SPP satellite, the 



precession rate increases by a factor of y Mj M to about —235 arcseconds per year. It is 
rather interesting that the SRP increases the rate of precession due to the oblateness of the 
sun, whereas we found that the precession rate due to the curvature of spacetime is decreased 
by the SRP. Perhaps the SRP could be used to disentangle these two types of precession. 

Since these two precessions occur in opposite directions, then for a certain choice of 
parameters they can actually cancel each other out. For a SSP satellite orbiting at 0.05 AU, 
for example, this would occur if the solar sail has a ~ 0.00163 kg/m 2 . Of course, it is more 
feasible to vary the average orbital radius for a fixed a until the precessions cancel. 
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2.5 Net electric charge on sun 

We will now consider the effect that a small amount of net charge Q on the sun would have 
on an SSP satellite with charge q. It has been suggested that the sun has a net charge of up 
to Q ~ 77 C [19]. The spacetime is then described by the Reissner-Nordstrom metric, which 
has the form ([H]), where the metric function / is now given by 

2GM Gk e Q 2 

f = 1 — + 4 2 , 56 

c r err* 

and the Coulomb constant k e = 8.988 x 10 9 N m 2 /C 2 . Taking into account the electric force, 
we find that 



2^3 



T 



2 



An 2 r 
GM 



k k e qQ k e Q 2 k 2 e q 2 Q 



1 + ^- + ^^ + -^ + 



(57) 

c 2 r GmM c 2 Mr G 2 m 2 M 2 _ 
where we have kept only terms up to order Q 2 . Note the Q 2 term in the period, due 
to the backreaction from the sun's charge on the geometry, is present even for a neutral 
SSP satellite. However, this term increases the period by an amount of only 10~ 35 s (10 -38 s 
without the SRP), which reflects the fact that the backreaction of the charge on the geometry 
is negligible, as to be expected. 

Primarily due to the photoelectric effect, but also the Compton effect and electron- 
positron pair production, a solar sail made from Beryllium, for example, will equilibrate to 
a charge per area of 0.065 C/m 2 [20,21]. For an SSP satellite with the specifications given 
in (jSj) and a mass of 1000 kg, this gives a charge g = 5x 10 4 C and an increase in period of 
about 230 s (0.05 s without the SRP). Thus, due to the SRP, even a small charge Q could 
certainly yield a measurable increase in the period, making this a potentially powerful test 
for net charge on the sun. 

2.6 Cosmological constant 

Since the SRP enhances a variety of small effects, including those associated with spacetime 
curvature, one could ask whether an SSP satellite could be used to test for the presence of 
a cosmological constant. In fact, supernovae observations suggest that our universe might 
have a very small positive cosmological constant [22,23]. In the presence of a cosmological 
constant A, the spacetime in the vicinity of the sun is described by a metric of the form (jUj), 
where the function / is now given by 

2GM Ar 2 
c 2 r 3 

Considering the effect of A for r ^> GM/c 2 , we get 

c 2 A 

GM ' 3GM 



f = 1 - — • (58) 



47r 2 r 3 ,_>.., 

I + ttt^ 3 ■ (59) 
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For A m 10~ 52 m~ 2 along with the specifications in (J5J), this leads to an insignificant increase 
in period of 10~ 17 s (10~ 21 s without the SRP). While one might be able specify values of 
a and r such that the correction factor is actually significant, the period would then be too 
large to make this a feasible test for the presence of a cosmological constant. 



3 Circular orbits out of the plane of the sun 

3.1 Review of non-Keplerian orbits 

We will now briefly review non-Keplerian orbits within the framework of Newtonian gravity 
[24-26,7,27]. The plane of a non-Keplerian orbit does not pass through the center of mass 
of the sun, and the SSP satellite is levitated above the sun, as shown in Figure EJ We work 
with spherical coordinates, where r is the heliocentric distance and 9 and <fi are the polar and 
azimuthal angles, respectively. We orient our coordinate system such that the SSP orbits 
around a circle at constant r and 9. Then the orbital radius is given by p = r sin 9. 











/ 



/ \ 

1 

reflected light 

Fi gure 2: Orienting a solar sail at an appropriate angle ip renders it possible to have a circular orbit outside 
of the plane of the sun. 

A non-Keplerian orbit can be maintained with a suitable pitch angle ip relative to the 
radially outgoing solar radiation hitting the surface of the sail. This enables one to control 
the direction of thrust due to the reflected portion of sunlight. By changing the pitch angle, 
one can also transfer the SSP satellite between different non-Keplerian orbits whose orbital 
planes have different orientations, as well as between non-Keplerian orbits and orbits within 
the plane of the sun. 

Note that while we are restricting ourselves to planar sails, for which the thrust is normal 
to the sail surface, the analysis for curved sails is similar as long as the sail curvature is 
symmetric about its center. The analysis below also carries over to non-Keplerian orbits of 
magnetic [28,29] and electric [30,31] sails, which rely on the solar wind rather than sunlight. 
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3.1.1 Perfectly reflecting solar sail 

For a perfectly reflecting solar sail (jj = 1), the magnitude of the acceleration due to the 
solar radiation pressure is given by 

KCOSIp 

a SR p = ■ (60) 

The factor of cosip is present because we have projected the SRP along the normal direction 
to the solar sail. We can break this acceleration into its components expressed in spherical 
coordinates as follows: 

a1 SRP = a SRP cos t/> , a e SRP = -a SRP sin ip , a% RP = , (61) 

where r is the heliocentric distance, 9 is the polar angle and <fi is the azimuthal angle. 
Including the contribution of acceleration due to the gravitational attraction of the sun gives 

K COS 2 1p — GM a K jl 

a = , a = — -coswsinw, a/ = . (62) 

Acceleration in terms of spherical coordinates can generally be written as 

a r = r — r6 — r sin 2 6 (f) 2 , 

a e = r9 + 2r8 — r sin 9 cos 9 <p 2 , 

= r - l d t (r 2 sin 2 00) , (63) 

where " = d t . Equating the components of acceleration in (1621 and (16"3"j) gives the equations 
of motion for the solar sail. For constant r and 9, these equations reduce to 

a = — r sin 9 <p = — , 

a 6 = —r sin 9 cos 9 <j) 2 = — - cos ip sin ip , 

a* = r - 1 d t (r 2 sm 2 9 <j)) = 0, (64) 

where 

M = M - ^ cos 2 ^ . (65) 

The (p equation implies that r 2 sin 2 9 <j) = Ln is a constant of motion, which corresponds to 
the angular momentum per unit mass of the solar sail. As before, we will use the subscript 
N to denote that this is a conserved quantity within the Newtonian approximation. 

We can write = 2tt/T, where T is the period of a complete orbit. Substituting this into 
the a r equation in (|64j) enables us to express the period in terms of the heliocentric distance 
and polar angle as 

T 2 = ii r 3 sin 2£ (66) 

GM 
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which bears some similarity to Kepler's third law, in that T ~ p 3//2 . From (1661) and the a 6 
equation in we find that 



GM 



+ 



An 2 GMr 3 cos 2 



tan ip r 



T 2 ' GMT 2 — 47r 2 r 3 sin 2 
2vr 2 r 3 sin 29 



GMT 2 — 4ir 2 r 3 sin 2 9 



(67) 



where the subscript r indicates that these are the parameters for a perfectly reflecting solar 
sail. Thus, the pitch angle ip r and the K r parameter are both determined by the period, 
heliocentric distance and polar angle of the orbit. For fixed r = 0.5 AU and T = 70 days, 
K r and ip r are shown in Figure [3]for < 9 < n/2. Note that for 7r/2 < 9 < ir, the plot of K r 
is flipped horizontally and that of ip r is flipped vertically. n r monotonically increases with 9 
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Figure 3: K r and tp r as functions of 6 for the SSP satellite. 

and reduces to the k given by (jl]) for 9 = n/2. i\) r is at a maximum for 9 = 7r/4 and vanishes 
for 9 = and rr/2. At first sight, it might seem odd that ip r decreases with 9 in the region 
7r/4 < 9 < 7r/2. However, from Figure [2] we see that the angle of the normal direction of the 
solar sail relative to the line passing through the poles of the sun is given by 9 — ip r , which 
monotonically increases with 9 down to the equatorial plane, as expected. 

3.1.2 Partially absorbing solar sail 

We will now consider the more realistic scenario in which a portion of the light is absorbed 
by the surface of the solar sail, which has also been considered quite a while ago in [24,25]. 
We saw above that reflected light pushes the solar sail at an angle ip relative to the radial 
direction. For a partially absorbing solar sail, the absorbed light pushes the solar sail radially 
outwards. The fraction of light reflected is 2t] — 1 and the fraction of light absorbed is 2(1— if). 
Putting all of the contributions to the acceleration together, we find that 

r (1 ~ V)"+ (2ry-l)« cos 2 jj-GM e (2r, - l)k . 
a = , a = cos ip sin ip , cr = , (68) 
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where k = L s /(2nca). This reduces to (|62|) for 77 = 1. Equating the components of acceler- 
ation in ( 1631) and ( 1681) and taking r and # to be constant gives 



a 



-r sin 2 9 d> 2 



GM* 



a = — r sin # cos 9 



(277 - l)k 



cos ^ sin ip . 



r 1 dt(r 2 sin 2 9 



where we now have the effective solar mass as 



M* = M-^ [(1 - r)) + (277 - 1) cos 2 if)] 



(69) 



(70) 



Plugging = 27r/T into the a r equation in ( 1691) yields the expression ( |66j) with M replaced 
by M*. Combining this with the a 9 equation in (|69|) gives 



1 - 



(2rj - 1)(GMT 2 - 4vr 2 r 3 sin 2 t 
4tt 2 (1 -r])r 3 sin(29) 

(GMT 2 - 4vr 2 r 3 sin 2 + tan 2 i/j a ) 
T 2 [77 H- (1 — rf) tan 2 




77(77 - l)(47r 2 r 3 sin 2^) 2 
(277 - 1) 2 (GMT 2 - 47r 2 r 3 sin 5 



(71) 



where the subscripts a indicate that a portion of the light is absorbed by the solar sails. 
Consider the case in which the solar sail is close to being a perfect reflector. Then e = 
1 — 77 1 and we have 



(tan -^a, n a ) = (tanipr, «v) x 



1 + 



(27r 2 r 3 sin 29)' 



(GMT 2 - An 2 r 3 sin 2 



(72) 



tanipr and K r are given by fIBTj) . As one might expect, the nonzero absorption of light means 
that we must have a larger pitch angle ip in order to have a non-Keplerian orbit. 

For fixed r = 0.5 AU, T = 70 days and r/ = 0.85, n a and ip a are shown in Figure H] for 
< 9 < 7r/2. The partial absorption of light by the surface of the solar sail causes K a to 
decrease and ip a to increase, which can be seen for 77 = 0.85 by comparing Figures [3] and HI We 
have confirmed that these changes in K a and ip a occur monotonically in 77 for 0.5 < 77 < 1. 
The 9 dependence of K a and ip a nevertheless remain qualitatively the same, regardless of 
whether any light is absorbed by the sail or not. This agrees with the conclusions of [24,25] 



3.2 Static curved spacetime 

We will now consider the effects of curved spacetime on non-Keplerian orbits. We will 
begin in this section by considering the sun to be a static central mass, and ignore the 
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Figure 4: K a and ip a as functions of 9 for the SSP satellite. 



frame dragging due to its rotation. Then the exterior spacetime of the sun is approximately 
described by the Schwarzschild metric, which is given by ([9]) and ( fTUl) . 

The acceleration 4-vector of an object moving within this background can be obtained 
by taking the covariant derivative of the velocity 4-vector = dx^/dr as follows: 



a 



dv^ 

"- — + r>V > r ^ = w + to* - g aP ,.) , (73) 



dr 



where x M = (t, r, 9, 0) and r is the proper time. We find that 



a r = f+^-/i 2 -^ T r 1 f 2 -/r(0 2 + sin 2 ^ 2 ) 



GM . 2 _ GM 
a 6 = rO + 29r — r sin 9 cos 9 <p , 

a* = r-^fr 2 sin 2 0) . (74) 

Note that the spatial components of the acceleration in (j74p reduce to those in (I63p in the 
Newtonian limit c 2 r >> GM. 

The acceleration due to the solar radiation pressure is given by 

(1 - rj)k + (277 - 1)kcos 2 ^ e (277 -1)k , • , * <a « 

a r = ^ -, a = - - 7 cos ib sin V> , a = a = . (75) 

Equating the components of the acceleration in ( 1741) and ( 1751) gives the equations of motion 
for the solar sail. For constant 9, these equations reduce to 

a' = f- l d T (fi) =0, 

a r = f + 9E. fP _ 9" V - /r si „> 9 = (l-))*+(fr-l)S<™V 



. „ . ; 2 (277 — 1)k 
a = — r sin 9 cos 9 <p = 

= r-^r 2 sin 2 # 0) = 0. (76) 



cos ip sin 
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The t and equations imply that 

ft = E/c 2 , r 2 sin 2 9 <p = L , (77) 

are constants of motion. The first integral of the r equation gives 

2 E 2 2 2GM* L 2 
r 2 = — -c 2 + -/, 78 

where M* is defined in ( 1701) . 

Circular orbits must satisfy the conditions given in ( ITBl) . which yields 

2 _„ 4 / 4GM-c 2 r \ c 2 GM* - c 2 GMV 2 
Using (1771) . we find the orbital period T to be given the expression 



2 



T 2 = ^Usin 2 #r 3 
GM* 



1 + k [(1 - 77) + (277 - 1) cos 2 V] 



c 2 r - 4GM) 



(c 2 r - 2GM) 2 

Keeping only the leading general relativistic correction gives 



An 



2 



sin 2 9 r 3 



1 + 



-r?) + (2/7- 1)cos 2 V] 



c 2 r 



(81) 



GM* 

From (|8ip and the a 9 equation in (170]) . we find the following equation for the pitch angle ip: 



(2r7 - l)c 2 tan ip(GMT 2 - 47r 2 r 3 sin 2 5) 
27r 2 r 3 sin 20[(1 - 77) (1 + tan 2 ip) + 2<q- 1] 



c 2 



3GM + 47r 2 r 2 sin 2 + 27r 2 r 2 sin 20[(1 - 77) (1 + tan 2 if)) + 2rj - 1] 
r T 2 (2i] — 1)T 2 tan if) 

Since all of the general relativistic corrections on the right-hand side go as c -2 , we can solve 
for if> perturbatively as if) = ip a + c~ 2 ip\. Then we find 

tanip = tan7/> a H — -(1 + tan 2 if> a ) tan^i, 
& 

f 2(2r] - 1) tan ^ a tan ^1 47T 2 r 2 sin 2 0\ 
K = Ka V 1+ c 2 [ V + (I-77) tan 2 ij a ] &T 2 J' (83) 



where tan if) a and n a are given in ( fflj) and 

27T 2 r 3 sin 20 [77 + (1 — 77) tan 2 ip a ] 



tan-0! 



(2/7 - 1)(1 + tan 2 ip a )(GMT 2 - 47r 2 r 3 sin 2 ^[77 - (1 - 77) tan 2 ifi a ] 

/GM 47r 2 r 2 sin 2 27r 2 r 2 sin 2% + (1 - 77) tan 2 \ 
X + T 2 + (277- l)T 2 tan^a J ' 
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For a perfectly reflecting solar sail, (153"]) reduces to 

47r 2 GMr 2 sin20 



t&nip = tan?/v + 



c 2 (GMT 2 — 47r 2 r 3 sin 



2 



/ GM(47T 2 r 3 sin20) 2 47r 2 r 2 sin 2 fl\ 

K ~ Kr [} + c 2 r[(GMT 2 ) 2 + 87r 2 r 3 sin 2 0(2vr 2 r 3 - GMT 2 )] c^T 2 J ' ^ 

where tan-?/v and k r are given in ( 1671) . We find that the curvature of spacetime has an 
extremely small effect on a and ip. For instance, for r = 0.5 AU, T = 70 days and 9 = 45°, 
a increases by about 4 x 10~ 13 kg/m 2 and tp increases by about 4 x 10~ 8 degrees. 



3.3 Frame dragging 

3.3.1 Orbits parallel to the equatorial plane of the sun 

As is the case for orbits within the plane of the sun, non-Keplerian orbits will also be affected 
by frame dragging due to the sun's rotation. We will first consider the effect of frame dragging 
on non-Keplerian orbits which are parallel to the equatorial plane of the sun. 

The exterior spacetime of the sun is approximately described by the large-distance limit 
of the Kerr metric, which is given by (1331) . From (1731) . we find that the components of the 
acceleration 4-vector of an object moving within this background are 

a 1 = f l d T (ft) + —<P d T 

a r = r + ^/i 2 -^rV 2 -M^ + sin 2 ^ 2 )-^/#sin 2 ^ 

AG J 

a = rO + 26 f — rsin#cos# 2 H — — -0tsin#cos# , 

c 2 r 2 

a* = r^d^r 2 sin 2 9 <P) + ^i(^^r -29 cos . (86) 

Note that the reduce to what is given by (1T4"]) in the limit J = 0. Equating the components 
of the acceleration in (|86|) and (!75j) gives the equations of motion for the solar sail. For 
constant 9, the equations for a r and a 9 reduce to 

(1 — 7])k + (2l] — 1)kcos 2 l[) 

AG J • ■ (2ri — l)k 
a 9 = — r sin 9 cos 9 <p 2 H — — — -<j)i sin 6 cos 8 = cos ip sin ip . (87) 

The first integral of the r equation gives 

. 2 _ -E 2 c2 + 2GM* L 2 AGJEL sin 
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where M* is defined in ( 1701) . 

Circular orbits must satisfy the conditions given in (1181) . which yields 



E 



c 4 + 



/ AGM 



c 2 GM*r 2 



c 2 r 



- c 2 r\ c 2 GM* 
3GMJ 



2c 2 GJsm6 GM*X 



c 2 r 



3GM (c 2 r - 3GM) 



r r(c 2 r - 3GM) 2 

6GJ sin 9 /—^ — — 

VGM*rX . 



where 



X = c 4 r 2 - c 2 G(3M + M*)r + AG 2 MM* . 
Using (1361) and (|89|) . we find the orbital period T to be given by 



(89) 



(90) 
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Air 2 
GM* 



1 + 



sin 2 9 r 3 1 + h [(1 - tj) + (2rj - 1) cos 2 
2GJsm9 / L E 



(c 2 r - AGM) 
{c 2 r - 2GM) 2 



c 2 r 



r 2 E + c 2 fL 



(91) 



As we did for the case of orbits within the plane of the sun, we have presented T 2 as a product 
of three factors. The first factor is the same as for non-Keplerian orbits in the Newtonian 
approximation for gravity. The second factor is due to the simultaneous effects of the SRP 
and the curvature of spacetime outside of a static central body. The third factor embodies 
the combined effects of the SRP and frame dragging due to the rotation of the sun. 
Keeping only the leading correction due to spacetime curvature gives 



T 



An 2 
GM* 



sin 9 r" 



1 + 



k [(1 -rj) + (277 - l)cos 2 ?/>] 



c 2 r 



1 + 



2 VG J sin^ 



(92) 



As was the case for orbits within the plane of the sun, the SRP increases the effect that frame 
dragging has on the orbital period, since it is the renormalized mass M that appears in the 
last factor. Although the solar sail parameter R and the pitch angle ip could be determined 
in terms of the other parameters from fl92|) and the a 9 equation in flHTj) . frame dragging has 
a negligible effect on k and ip. 



3.3.2 Polar orbits outside of the plane of the sun 

We will now consider the effect of frame dragging on non-Keplerian orbits which are parallel 
to polar orbits, and outside of the plane of the sun. The plane of non-Keplerian polar orbits 
undergoes the Lense-Thirring effect, as shown in Figure [51 

We do not present the resulting equations since they are rather cumbersome. However, 
recall from the discussion of frame dragging for polar orbits within the plane of the sun that, 
up to linear order in J, the precession frequency was not affected by the SRP nor the polar 
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Figure 5: The Lense-Thirring effect for non-Keplerian polar orbits. 



orbital motion of the satellite, and depended only on the heliocentric distance r and not the 
polar angle 9. Up to leading order in J, these characteristics also hold for non-Keplerian 
orbits and the precession frequency is given approximately by (j46j) . Thus, in one orbital 
period given approximately by flHTj) . the angle of precession is 

. , 47rGJsin0 . , 

A6 w . (93) 

c 2 VGM*r 3 

3.4 Oblateness of the sun 

We will now consider the effect of the quadrupole mass moment of the sun on non-Keplerian 
orbits. The gravitational potential is given by f j4"8l) . from which we find that the components 
of acceleration due to the gravitational potential are 

GM 3GMJ 2 R 2 . 2/1 3GMJ 2 R 2 , nn . , 
n r = (3cos 2 #-l), a e = _J_siii(20), a* = . (94) 



r 
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Including the contribution of acceleration due to the SRP gives 

GM* 3GMJ 2 R 2 



r 2 2r 4 



(3 cos 2 9 



a 6 = -U—J-anW) ^— sin(20) , 

a? = 0, (95) 

where k = L s /(2ttc<7) and M* is given by flTUl) . Equating the above components of accelera- 
tion to those in f[6"3l and taking r and 9 to be constant gives 

• 2, 12 GM* SGMJ 2 R 2 2 
a = —r sin 9 6= (3 cos 9 — 1) , 

r 2 2r 4 

o ■ a a a? ^ ~ ^ k ■ 3GMJ 2 R 2 . 
a = — rsm9cos9 <p = sm[2ip) sm(20) , 

-- r- 1 d t {r 2 sin 2 9 <p) = 0, (96) 



a' 



Substituting 6 = 2n/T into the a r equation yields 



„ 4vr 2 r 3 sin 2 9 f 3MJ 2 R 2 ,„.,,„, 

T = 1 + — =-= — l-3cos 2 # . 97 

GM* V 2M*r 2 1 
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Note that this reduces to the corresponding terms in (1501) when 9 = tt/2. Interestingly 
enough, the oblateness of the sun increases the period for orbits within the angular range of 
approximately 55° < 9 < 90°, whereas the period is decreased for 9 < 55°. 
From the a r and a 6 equations in fl96|) . we find 



,n 2 n 47r 2 r 3 sin 2 fl 3GMJ 2 R 2 ,„ 2n „ N 
k[1 - + (2^-1) cos 2 ^] = GM — + -^—(3 cos 2 0-1), 

4vr 2 r 3 3GMJ 2 R 2 * 



(2 ?7 -l)Ksin(2V) = 1-^ -^Jsin(2£), (98) 

from which one can express «; and the pitch angle ip in terms of the orbital parameters and 
77. Since these expressions are rather long, we present the idealized case in which the solar 
sail is a perfect reflector: 

3GMJ 2 R 2 r 2/1 A 2 , 2 , „, /GMT 2 + 7r 2 r 3 (7cos 2 #-5) 

1 -3cos 2 # + 47rVsin 2 (2#) ' 



2r 2 



(GMT 2 - 47r 2 r 3 sin 2 



tan^ ~ taxiipr 



3GMJ 2 R 2 T 2 (GMT 2 + 2vr 2 r 3 (5 cos 2 9 - 3)\ 
47r 2 r 5 V GMT 2 — Air 2 r 3 sin 2 J 



(99) 



where we have included the leading term in J 2 <C 1. We find that the oblateness of the sun 
has an extremely small effect on a and if}. For instance, for r = 0.5 AU, T = 70 days and 
9 = 45°, a decreases by about 2 x 10~ 12 kg/m 2 and ^ decreases by about 2 x 10~ 7 degrees. 



4 Conclusions 

We have considered the effects of various phenomena along with the solar radiation pressure 
(SRP) on bound orbits of a solar sail propelled (SSP) satellite. The consideration of the 
SRP on its own leads to a renormalization of the effective solar mass. However, when the 
SRP is coupled with other effects, such as spacetime curvature, then the resulting orbital 
characteristics depend on the solar mass M and the solar sail parameter k separately, and 
not simply in the combination of a renormalized mass. 

For simplicity, we first considered the effects on the period of a circular orbit in the plane 
of the sun. In the table below, we summarize our results for the change in period AT of an 
SSP satellite with the specifications given in (jSJ). 
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Phenomenon 


AT without SRP 


AT with SRP 


Static Curvature of Spacetime 


s 


0.6 s 


Frame Dragging 


4 x 10" 5 s 


0.01 s 


Quadrupole Mass Moment of Sun 


0.02 s 


105 s 


Octupole Mass Moment of Sun 


10" 7 s 


6 x 10~ 4 s 


Net Charge of Sun 


0.05 s 


230 s 


Cosmological Const. 


lO" 21 s 


10" 17 s 



Recall that the base period without the SRP is about 4 days, and with the SRP is about 
70 days. The SRP generally augments the change in period due to various phenomena by 
a factor of about 1000 or more. This might be sufficient to observe the AT from some 
phenomena, such as the quadrupole mass moment of the sun and a possible net charge on 
the sun. 

Although we have isolated the effects of each of these phenomena, the reality of course is 
that they all have simultaneous effects on the period. For instance, to compute the combined 
effects of frame dragging and the mass multipole moments of the sun on orbits, one could use 
the metric which describes a rotating deformed mass with all mass multipole moments [32]. 
Naturally, the challenge is to isolate the various effects from one another, especially since 
the SRP would also enhance additional effects that we have not discussed, such as the 
gravitational effects of the planets. 

As an example, consider the effect on the period due to the static curvature of spacetime 
and the SRP. Note that this has no counterpart without the SRP. Namely, when the SRP is 
not taken into account, the static curvature of spacetime does not alter the period as given by 
Newtonian gravity. Unfortunately, this rather interesting effect could likely be overshadowed 
by the dominant effect of the quadrupole mass moment of the sun. However, the effect of 
the mass moments of the sun depends on the plane of orbit, whereas the effect of the static 
curvature of spacetime does not. Thus, considering different orbital planes could aid in 
isolating the weaker effect of spacetime curvature. We have also confirmed that the higher 
mass moments of the sun have negligible effect compared to that of the static curvature of 
spacetime. 

Likewise, for a charged satellite, the effects of net charge on the sun could be smeared out 
by the sun's magnetic field. This too is unfortunate, since the magnitude of the net charge 
on the sun is not known very accurately, so any improved measurements of this are highly 
desirable. Since the magnetic force depends on the velocity of the satellite, the situation 
might be substantially improved by considering orbits with various orientations and speeds. 
However, since the magnetic field of the sun is highly variable, extracting the effects of net 
electric charge remains a universal problem. 
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We also considered orbits which deviated from being perfectly circular. There are two 
sources of perihelion shift. The curvature of spacetime causes the perihelion to be shifted 
in the direction of the orbit, whereas the oblateness of the sun causes the perihelion to be 
shifted in the reverse direction. Without the effects of the SRP, the perihelion shift due to 
spacetime curvature is generally larger than that due to the oblateness of the sun. However, 
this is not the case when the effects of the SRP are taken into account on an SSP satellite. As 
shown in the table below, the SRP decreases the perihelion shift due to spacetime curvature 
but increases the perihelion shift due to the oblateness of the sun. 



Cause of Perihelion Shift 


A0 without SRP 


A0 with SRP 


Spacetime Curvature 


707year 


4'Yyear 


Oblateness of Sun 


-14" /year 


-235'Vyear 



Thus, SSP satellites might enable one to disentangle the perihelion shift due to the oblateness 
of the sun from the general relativistic perihelion shift by making the former one more 
pronounced. Of course, these are both small effects which would need to be isolated from a 
plethora of other types of effects. 

We have also investigated circular orbits outside of the plane of the sun, which are 
commonly referred to as non-Keplerian orbits. Such orbits are possible for a solar sail with 
an appropriate design parameter k and pitch angle ip, both of which depend on the specified 
orbital period, radius and distance from the plane of the sun. We have computed how k 
and ijj depend on the heliocentric distance and polar angle of the orbit, as well as on partial 
absorption of light by the surface of the solar sail. We have also considered how k and ip are 
altered due to general relativistic effects as well as the oblateness of the sun. 

For both orbits within the plane of the sun as well as non-Keplerian orbits, we have 
considered the Lense-Thirring effect, which is the precession of polar orbits. To leading 
order, the rate of precession is not altered by the SRP. Namely, any satellite in polar orbit 
at a distance from the sun of 0.05 AU will experience a rate of precession of about 0.03 
arcseconds per year. 

For the case of perfectly reflecting solar sails, types of non-Keplerian orbits have been 
classified and their linear stability characteristics have been studied [7] . In particular, stable 
and unstable regions in space have been plotted. It would be interesting to see how the 
boundary between stability and instability is altered due to various effects such as the partial 
absorption of light, the oblateness of the sun and curved spacetime. The patching between 
two non-Keplerian orbits of different orientations, as well as between a non-Keplerian orbit 
and a Keplerian orbit, has also been studied [7]. It would be interesting to see if the allowed 
types of patching are changed by these various effects. 
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In this study, we have assumed that the intensity of solar radiation goes as the inverse 
square of the heliocentric distance. However, this assumption is not valid when the finite 
angular size of the sun is taken into account- along with the oblateness of the sun, sunspots, 
etc. Although the deviation from an inverse square law is certainly small, orbits close to the 
sun as well as long-duration orbits will experience a cumulative perturbation. 
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